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ABSTRACT 


Global  temperatures  are  obtained  using  a  linear  least 
squares  regression  method  with  satellite  radiation  measure- 
ments, particularly  those  obtained  from  the  SIRS-B  aboard 
the  NIMBUS  IV  satellite.   Regression  equations  relating 
temperature  to  spectral  radiance  observations  are  employed 
The  regression  equations  were  then  applied  to  independent 
observations  to  determine  the  feasibility  of  temperature 
determinations  over  sparse  data  regions. 
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X.   INTRODUCTION 

The  launch  of  the  NIMBUS  IV  on  8  April  1970  was  the 
second  step  in  a  process  of  developing  a  means  of  determin- 
ing the  three-dimensional  state  of  the  atmosphere  by  indir- 
ect means.   The  orbit  of  the  NIMBUS  IV  was  a  near-circular, 
polar-orbit,  with  a  height  of  approximately  600  nm .   The 
new  instrument,  known  as  SIRS-B  (which  is  similar  to  the 
instrument  on  NIMBUS  III,  known  as  SIRS-A))  has  the  added 
capabilities  of  "side  looking"  scan  for  greater  area  cover- 
age, and  the  determination  of  the  tropospheric  water  vapor 
content . 

The  calibration  procedure  used  to  transform  the  radia- 
tion sensed  in  the  eight  channels  of  the  SIRS-B  to  spectral 
radiances  is  discussed  in  the  NIMBUS  IV  USERS  GUIDE  (1970). 

The  temperature  profile  of  the  atmosphere  within  the  field 

2 
of  view,  comprising  an  area  of  (225km)   at  the  subsatellite 

point,  is  deduced  using  the  channel  spectral-radiance  en- 
semble.  The  particular  wavenumbers  sensed  are  small  inter- 
vals in  the  atmospheric  window  and  the  15-micron  band  of 
C0_.   The  channel  numbers  and  wavenumbers  are  listed: 

CH.  .NO.     1       2  3      4      5      (3      7      8 

l 

WAVE.  NO.    899.0  750.0  734.0  709.0  701.0  692.0  679.8  668.7 
v  (cm"1) 


Through  the  equation  of  radiative  transfer,  the  meas- 
ured radiances  are  related  to  the  temperature-structure 
T(p)  of  the  atmosphere,  the  temperature  at  the  earth's  sur- 
face and  that  at  the  cloud  tops  which  may  be  within  the 
f ield-of-view .   Assuming  that  a  single  cloud  layer  exists, 
the  following  set  of  integral  equations  can  be  constructed 
for  the  measured  radiance  N(v.)  in  channel  i: 

T 

N(v.)=A  j  B[v.,T(pc)]T(vi,pc)  -  /  CB[vi,T(p)]dT(vi,p)  ) 

(1) 

s 

+  (1-A)  (  B[v.T(ps)]  T(v.,ps)  -J       B[vi,T(p)]dT(v.,p)  ) 

i=l,  2,  .  .  ,  ,  8  is  the  channel  number. 

N(v.)  =  spectral  radiance  at  wavenumber  v.. 
i  i 

B[v.,T(p)]  =  Planck  radiance  at  wavenumber  V   and  tempera- 
ture, (Eq.  (3)  below). 
i(v.  jp)3  fractional  transmi 1 1 ance  of  the  atmosphere 
in  the  spectral  interval  centered  at  wave- 
number  v.  and  assumed  known  as  a  function  of 
l 

the  pressure  p  at  the  radiating  level. 

A  =  the  product  of  the  fraction  of  cloud  cover 

within  the  f ield-of-view  and  the  cloud  emis- 

sivity  (which  is  usually  assumed  to  be  unity) 

Subscripts  c  and  s  refer  to  cloud-top  and  earth  surface. 

The  special  case  of  equation  (1)  where  A=0  gives  the 

clear-column  or  corrected  radiance  N  (v.).   This  term  is 

c   1 

contained  within  the  second  brace  of  equation  (1).   It  is 
this  quantity  which  must  be  determined  in  all  channels  in 
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order  to  convert  the  corrected  radiance  set  into  a  vertical 
temperature  profile  T(p).   The  clear-column  spectral  radi- 
ance at  the  top  of  the  atmosphere  at  wavenumber  V.  is  there- 
fore : 

Nc(v.)  -  B[vi,T(ps)]T(v.,ps) 

Ts  (2) 

-  /  B[vi,T(p)]  dT(v.,Ps) 

where 

-1 


B 


[vifT(p)]  =  c±V3±    jexp  [c2V./T(p)j  -  lj         (3) 

is  the  Planck  function  at  wavenumber  v.  . 

l 

Here  T=T(p)  is  the  temperature  (assumed  to  be  in  thermo- 
dynamic equilibrium  at  level  p)  ,  and  c ..  and  c  _  are  known 

cons  tan  ts 

-5        2     -1      -1 
c   =  1.9061  x  10    erg  cm   sec    ster 

c2  =  1.43868  cm  °K 

T(v.,p)  is  the  transmi t tance  at  wavenumber  V.  of 
l  l 

the  atmosphere  above  level  p. 

Smith,  Woolf  and  Jacobs  (1970)  have  used  equations  (2) 

and  (3)  to  de  f  ine  a  weighted-mean  black-body  temperature 

T„  .  in  terms  of  the  corrected  radiance  N  (v . )  which  would 
B ,  l  c   i 

be  observed  from  an  atmosphere  of  structure  T=T(p) 

Nc(v.)  =  B(V.,TB.)  W 

Here  T    is  independent  of  pressure.   Solving  the  relation- 
B  i 

ships  (3)  and  (4)  for  T„ .  gives 

B  i 

T  (v  )  =  c?  v  /In  (  [c  v.3/N  (V,)]  +  1  }  (5) 

Smith  et  al   (1970)  have  provided  multiple  regression  equa- 
tions which  relate  NMC  effective  grid  point  temperature 
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T(p.)  at  SIRS-A  scan  spots  to  the  sixteen  predictors  T  (v . ) 
J  B   i 

2 
and  T  (v . )  .   They  have  done  this  by  relating  a  two-week 

NMC  data  file  of  temperatures  at  standard  pressure  levels 

p.  to  the  corrected  radiances,  and  arrived  at  an  equation 


of  the  form 


S 


T(p.)  =  T(p  )  +  I       a(v,,pJ[T,(v,)-  Tjv,)] 
J        J    i  =  l 


S 


i,rj/l  Bx  f    B  v  i 

2 


+  Z   a'(v.,p  )  [TB(V.)-  TB(V.)] 
1  =  1         J 


(6) 


A  similar  study  will  be  made  for  T(500)  and  T(300)  using 
a  data-base  of  corrected  radiances  from  SIRS-B  and  NESC  de- 
rived temperature  profiles  T(p.)  together  with  simultaneous 
FNWC  ef  f e  c  t ive  grid-point  soundings  (available  from  the 
Fleet  Numerical  Weather  Central  ADP  tapes)  at  locations  of 
the  scan-spots  in  both  space  and  time.   For  the  SIRS-B  data 
file,  a  revised  NESC  formula  similar  to  equation  (6)  was 
available  for  the  conversion  of  corrected  radiance-sets 
into  a  statistically  consistent  profile  T=T(p). 
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II.   DATA  PROCESSING 

A.   THE  ORIGINAL  DATA 

The  original  data  utilized  in  this  paper  were  provided 
by  the  National  Environmental  Satellite  Center  (NESC) ,  and 
consisted  of  computer  listings  of  a  selected  set  of  SIRS-B 
scan  spots  recorded  during  December  1970.   In  addition  to 
the  geographical  coordinates  and  time  of  each  observation, 
there  were  also  listed  the  raw  (uncorrected)  radiances  in 
each  channel.   The  raw  radiances  are  assumed  to  be  derived 
from  equation  (1)  for  the  one  cloud-layer  model  and  the 
corrected  radiances  resulted  from  the  systematic  elimina- 
tion of  cloud  contamination  effects  from  each  channel  [fol- 
lowing a  technique  of  Smith   et  al   (1970)]. 

The  atmospheric  window  channel  (channel  1)  was  inopera- 
tive during  this  period;  therefore,  the  values  of  the  cor- 
rected window  channel  were  chosen  from  locations  where  the 
surface  temperature,  based  on  the  NMC  data  file,  was  known 
to  be  at  a  pressure-level  close  to  p  =  1000  mb ,  where  T   is 
known.   Furthermore,  since  the  radiance-correction  computa- 
tion involves  iterative  use  of  equations  (2)  and  (6) (using 
an  equivalent  "clear-column"  temperature  profile,  T=T(p), 
at  each  iteration),  all  original  radiance-sets  chosen  could 
thus  be  "tied  down"  by  the  selection  of  scan-spots  having 

surface  pressure  values  p   =  1000  mb ,  with  T   known. 

r  s  s 
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It  should  be  noted  that  the  data  samples  provided  were 
enhanced  by  the  side-viewing  capability  of  the  SIRS-B.   The 
side-viewing  capability  enabled  the  scanning  of  more  atmos- 
phere-earth columns  both  on  and  off  the  subtrack  for  which 
the  approximation  p   =  1000  mb  was  valid.   Slant-columns 
thus  viewed  are  useful  for  reduction  to  clear-column  tem- 
perature profiles  T  =  T(p)  after  the  slant  path  is  normalized 
to  a  single  atmospheric  depth  over  the  spot  viewed. 

In  general,  virtually  all  scan  spots  have  some  effective 

2 
cloud  cover  (A  >  0)  in  view  of  the  large  area  (225  km)   be- 
ing viewed.   Of  the  seven  C0„  channels  being  used,  the  two 
channels  subject  to  the  maximum  attenuation  by  the  inter- 
position of  a  "cool"  cloud  top  within  the  desired  T  =  T(p) 
profile  are  channels  2  and  3. 

Smith,  Woolf,  and  Jacobs  (1970)  have  indicated  the  pro- 
cedure involved  for  correcting  for  cloud-contamination  of 
the  eight  channels  of  radiances.   Based  upon  the  actual  NMC 
climatological  update  of  each  of  the  standard  pressure- 
levels,  (10,  30,  50,  100,  150,  200,  250,  300,  400,  500, 
700,  850,  and  1000  mb)  they  have  derived  multiple  regression 
equations  relating  T(p.)  to  the  clear-column-radiance 
N  (v.),  or  to  the  blackbody  temperature-equivalent  of  the 
clear-column  radiance.   For  example,  Smith  et  al   (1970) 
obtained  regression  equations  essentially  of  the  form  of 
equation  (6)  and  have  listed  the  least-squares  best-fit 
coefficients  obtained  for  a  period  in  September  1969  appli- 
cable in  a  zonal  band  35-55  N.   As  a  second  approximation, 
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Smith  et  al   have  also  generated  a  six-predictor  analog  of 
equation  (6)  with  channels  2  and  3  deleted  for  the  same 
time  frame  as  utilized  for  the  eight-predictor  equation. 
In  both  cases,  they  have  deduced  equation  (6)  and  its  six- 
predictor  analog  by  employing  a  large  sample-size  of  simul- 
taneous values  of  corrected  radiance-sets  and  of  values  of 
T(p.)  (at  standard  pressure  levels),  interpolated  from  the 
NMC  file  to  the  location  of  the  SIRS  scan  spot.   The  six- 
predictor  and  eight-predictor  regression  equations  are  con- 
tinually updated  to  a  revised  two-week  data  base  after 
three  additional  data-days  are  accumulated.   Such  equations 
were  also  available  during  the  SIRS-B  period  of  December 
1970. 

The  procedure  for  correction  of  the  i_th  channel  radi- 
ance makes  iterative  use  of  the  updated  SIRS-B  operational 
six-predictor  regression  equation  in  deducing  a  first  guess 
of  T=T(p).   At  the  first  guess,  the  uppermost  five  radiances 
(N,,...,N  )  which  are  least  affected  by  cloud  contamination 

H  O 

are  used  as  the  first  guess    correct-radiance  set,  together 
with  the  prescribed  N-  from  the  window-channel  blackbody 
"fix"  at  the  scan  spot.   With  these  six-predictors  known,  a 
first  guess  of  the  temperature  profile,  T(p.),  at  all 
standard  levels  is  computed.   After  the  first-guess  of  the 
temperature  profile  is  made,  the  profile  is  then  applied  to 
channels  2  and  3  to  determine  the  "best  fit"  for  cloud 
height  and  effective  cloud  cover  A  using  equation  (1),  re- 
written for  simplicity  as 
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N(V  )  =  (1-A)  N  [V   T(p),p  1  +  AY  [V.,T(p),p] 

X  C     X  o  X 


(7) 


Here  Y  represents  the  contents  of  the  first  brace  of  equa- 


tion (1).   Also  c  (v . )  is  defined  as 
c(vi)  =  Nc  (v±)    -  NCVj.) 


(8) 


for  the  correction  to  be  applied  to  the  radiances  for  chan- 
nels 4  through  8.   Writing  equation  (7)  for  channels  2  and 
3  in  terms  of  the  corrected  values  of  N   based  on  T(p)  (the 
clear-column  temperature)  and  the  corresponding  computation 
of  Y  from  the  top  of  the  clouds  gives  two  simultaneous 
equations.   From  these  two  simultaneous  equations,  one 
solves  for  the  two  unknowns,  A  and  p  ,  using  indirect 
methods  for  solving  for  p   (i.e.  making  separate  tests  for 
p   =  200,  250,  ...,  850  mb)  and,  in  turn,  determining  the 
value  of  A  which  "best-fits"  equation  (7)  expressed  in 
terms  of  N(v»)  and  N(v_).   These  values  of  A  and  p   are 
substituted  into  equation  (7)  for  a  second  corrected  esti- 
mate of  N ,,.... ,N_and  a  new  temperature  profile  is  calcu- 
lated.  The  iterative  procedure  is  continued  in  this  manner 
and  usually  converges  within  four  iterations,  i.e.,  the 
correction  c  (v  .  )  of  (8)  usually  converges 

c  (v  .  )  -*■  c  (v  .  ,  limit) 
i        i 

for  a  suitable  tolerance  within  four  iterations.  At  this 
point,  the  finalized  temperature  profile  T(p)  may  be  used 
to  compute  the  clear-column  radiances  for  channels  2  and  3 
as  well  as  the  five  uppermost  CO-  channels  using  Eq .  (1). 
Cloud  influence  is  thus  compensated  when  the  final  choice 
of  temperature  profile  has  been  made. 
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In  a  sizeable  fraction  of  the  cases  selected,  the  con- 
vergence c(v.)  required  more  than  four  iterations  which 
means  that  a  relaxed  requirement  for  convergence  had  to  be 
employed.   This  generally  occurred  with  cloud  tops  at  or 
above  500  mb .   Such  cases  of  below-normal  attenuation  re- 
covery were  indicated  on  the  NESC  data-listings  as  decreased 
reliability.   Nevertheless,  no  data  cases  were  discarded 
for  such  reasons. 

B.   ADJOINING  OF  FNWC  DATA 

Analyzed  temperature-field  data  for  500  and  300  mb  cover- 
ing the  period  from  00  GMT  24  December  1970  through  00  GMT 
29  December  1970  was  provided  by  Fleet  Numerical  Weather 
Central  (FNWC) .   The  FNWC  temperature  data  had  not  been 
correlated  with  radiances  as  was  the  NESC  data.   The  FNWC 

temperature  fields  at  500  and  300  mb  [T_(5)  and  T  (3)]  were 

r  r 

spatially  and  temporally  interpolated  to  obtain  consistency 
with  the  SIRS  soundings.   Spatial  interpolation  was  accom- 
plished using  Bessel's  central-difference  interpolation 
[Haltiner,  1971]  to  the  latitude  and  longitude  of  the  SIRS-B 
soundings.   Temporal  interpolation  to  the  time  of  the  scan- 
spot  was  accomplished  utilizing  linear  interpolation  between 
temperature  fields  preceding  and  succeeding  the  time  of  the 
SIRS  sounding.   These  methods  were  employed  to  obtain 

statistical  dependence  of  1^(5)  and  T_(3)  upon  the  corrected 

r  r 

radiance  set  N. N„.   This  necessitated  the  printout  of 

1  o 
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FNWC  interpolated  temperatures  in  a  format  consistent  with 
the  requirements  of  the  program  BIMED  02R,  to  be  used  in 
(III). 

C.   DATA  DESIGNATION 

The  original  NESC  data  comprised  some  207  to  329  scan- 
spots  per  day  and  were  identified  by  the  specific  time  of 
the  SIRS  soundings  and  their  latitude  and  longitude.   The 
original  FNWC  data  tapes  consisted  of  temperatures  at  the 
levels  300  and  500  mb ,  at  each  grid  point  of  the  63  x  63 
Northern  Hemisphere  FNWC  grid.   The  FNWC  data  at  scan  spots 
were  obtained  by  interpolation  as  previously  described. 

A  three-day  data  base  of  both  NESC  and  FNWC  data  was 
thencompiled  using  data  from  0600  GMT  to  1800  GMT  for  24, 
25,  and  26  December  1970.   A  four-day  data  base  was  com- 
piled by  the  addition  of  data  from  0600  to  1800  GMT  27 
December  1970  to  the  three-day  data  base.   Data  sets  for 
independent  regression  tests  were  compiled  using  data  cen- 
tered on  00  and  12  GMT  27  December  and  00  and  12  GMT,  28 
December  (i.e.,  data-base  times  plus  12  and  24  hours). 

All  data  sets  were  divided  into  three  bands.   These 
bands  were  bounded  by  latitudes:  20-40  N,  band  1;  35-55  N, 
band  2;  and  50-70  N,  band  3.   The  bands  were  overlapped  to 
reduce  the  possibilities  of  discontinuities  across  the 
boundaries  when  regression  calculations  were  made. 
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The  designation  of  dependent  and  independent  variables 
was  made  and  the  dependent  variable  was  designated  T(i,j) 
and  the  independent  variables  were  designated  as  shown: 


CH.  NO.   12345678 

WAVE  NO.  899.0  750.0  734.0  709.0  701.0  692.0  679.8  668.7 

VAR.  NO.  N.      N0     N.     N.      Nc     N,     N_     N0 
12345678 

For  each  scan  spot  selected  in  the  test  period,  the 

time  and  location  were  encoded  on  IBM  cards,  as  were  the 

temperatures  T  (5)>  T.T(3),  T_(5),  T_,(3).   Here,  the  sub- 
N         N         r         r 

script  N  indicates  an  NESC-generat ed  temperature,  using  the 
convergent  procedure  described  previously,  and  the  subscript 
F  denotes  interpolated  FNWC  temperatures.   Also  encoded 
onto  IBM  cards  were  the  eight  channels  of  corrected  radi- 
ances (  provided  by  NESC  using  the  method  described  in  II. A.) 
for  each  of  the  scan  spots.   All  data  were  encoded  in  a 
manner  consistent  with  the  requirements  of  the  stepwise 
regression  program,  BIMED  02R  (Dixon  1966)  which  will  be 
discussed  later. 


19 


III.   THE  REGRESSION  METHOD 


A.   THE  STEPWISE  REGRESSION  PROCEDURE 

Many  problems  in  research  require  the  extensive  analysis 
of  large  amounts  of  data.   The  data  handling  process  should 
be  made  as  automatic  and  rapid  as  possible.   The  appropriate 
tools  for  analyzing  large  data  samples  with  numerous  inde- 
pendent variables  often  require  multivariate  statistical  re- 
gression techniques  to  attack  the  problem.   The  Biomedical 
Computer  Programs  (Dixon  1966)  were  developed  to  provide 
the  techniques  necessary  for  the  statistical  and  mathemati- 
cal analyses  required.   One  of  these  programs  is  the  Step- 
wise Regression  Analysis  Program  BIMED  02R. 

BIMED  02R  computes,  in  a  stepwise  manner,  a  sequence  of 
multiple  linear  regression  equations.   One  variable  is  added 
to  the  regression  equation  at  each  step.   The  variable  added 
is  that  which  accomplishes  the  greatest  reduction  in  the 
previously  unexplained  sum  of  squares.   This  variable  is 
also  the  one  which  has  the  highest  partial  correlation  with 
the  dependent  variable  at  the  particular  step  in  the  analy- 
sis of  the  variance.   Equivalently  it  is  the  variable  which 
would  have  the  highest  F-value  if  it  were  added.   The  F- 
value,  or  F-statistic  upon  entry,  F,  is  expressed  at  step  k 
as  (Dixon  1966) 


Fk(l,n-k-l)  - 


%(C.E.V.k)  -  %(C.E.k-l) 
%(U.E.V.  ,k) 
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where 

%(C.E.V.k)  is  the  percent  cumulative  explained  variance, 

step  k. 

%(C.E.V.k-l)  is  the  percent  cumulative  explained 

variance,  step  k-1. 

%(U.E.V.k)  is  the  percent  unexplained  variance  remain- 
ing at  s  tep  k . 

This  study  utilized  a  statistical  model  expressible  in 

the  form 

T    =    C    +C.N1+C0N_  +  C0N0  +  C/N.+C,Nc  +  C,N,+C^N.7  +  CoNQ  (9) 

o       11       22       33       4    4       55       66       77       88 

where  C   through  C n  are  to  be  determined  by  the  stepwise 
o       °    8  r 

least  squares  technique.   Equation  (9)  may  be  expressed  in 

the  more  compact  form  utilizing  the  matrix-product  equation 

l 
Ni 

t  =  (c  ,  c. ,  ,  cR)(:  )  (io) 

01  8   N 

W8 

In  (10),  the  notation  T  (3,5)  will  be  used  when  the  NESC 
data  at  300  or  500  mb  (respectively)  are  correlated  against 

the  predictors  (N1  , >NQ);  whereas  T_,(3,5)  is  the  nota- 

tion  for  the  corresponding  predictands  with  FNWC  tempera- 
tures.  Equation  (10)  may  also  be  used  in  application  with 
independent  data  occurring  up  to  12  or  24  hours  later  than 
the  dependent  sample  which  yielded  the  particular  coeffi- 
cient matrix  (C  ,  C.  ....,  CQ). 

o    1         o 

B.   RELATED  STATISTICAL  PARAMETERS 

Available  as  outputs  from  the  BIMED  02R  program  are 
'several  related  statistical  parameters.   Among  these  sta- 
tistical parameters  are: 

21 


a.  multiple  R 

b.  standard  error  of  estimate,  S.  E. 

c.  mean  value,  T 

d.  standard  deviation  a 

e .  F-value 
2 


f  .   R 


2  . 


The  significance  of  R   in  relationship  to  S.  E.  is  as 
follows  (Crow  et  al   1955) : 

(S.E.)2  -  a2  [  (n-l)/(n-k-l)]  (1-R2) 
where 


(11) 


n        _  9 
O2  =  I        (T.  -  T)     /n 
i  =  l    1 

is  the  variance  of  the  temperature  sample,  and 

n  =  sample  size 

i  =  sample-element  index 

k.  =  number  of  predictors 
Furthermore,  R  is  the  multiple  correlation  coefficient  after 
k  predictors  are  added,  with  emphasis  on  k=8,  in  order  to 

derive  maximum  information  from  the  S IRS-sounding .   It  fol- 

2 
lows  that  the  percentage  unexplained  variance  (1-R  )  for 

sample  sizes  in  the  range  171  to  228  as  contained  in  the 

dependent  data  sample  are  very  closely  approximated  by 

(S.E. /a)2  =  1-R2 

2 
The  percentage  explained  variance  is  then  given  by  R  ;  that 


is 


R2  =  1  -  (S.  E./a)2 


(12) 


since  (n-k-1) / (n-1 )  is  close  to  unity  for  the  dependent 
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sample  sizes  considered  here.   The  factor  (n-k-1) / (n-1) 
becomes  (n-2)/(n-l)  =  1  for  the  independent  test  cases 
considered  since  in  this  case,  there  is  only  one  way  of 
forming  the  predictor  to  be  listed. 

C.   THE  REGRESSION  PROCEDURE  APPLIED  TO  THE  DEPENDENT 
SAMPLES 

The  test  on  the  dependent  data  was  made  using  the  least- 
squares,  stepwise  regression  technique  of  BMD  02R  to  develop 

the  coefficients  C  ,  C.  ,....,  C0  .   The  coefficients  for  the 

o   1         8 

three-day  and  four-day  dependent  samples  for  the  T(5)  and 
T(3)  regression  equations  were  computed  by  use  of  this 
technique.   This  test  was  performed  on  both  the  NESC-regres- 
sion  generated  data  and  the  interpolated  FNWC  data.   The 
coefficients  for  the  NESC  three-day  data  sample  are  listed 
in  Table  1,  part  (a)  for  500  mb ;  and  part  (b)  for  300  mb . 
Table  2  contains  similar  results  for  the  four-day  sample 
base  of  NESC  data.   The  coefficients  for  the  FNWC  data  are 
listed  in  Table  3  (three-day  case)  and  Table  A(four-day). 

The  relevant  statistics  for  the  tests  on  the  several 
dependent  data  bases  are  listed  in  Table  5  and  6.   The 

statistics  included  in  these  tables  are  the  mean,  standard 

2 
deviation,  R  ,  and  standard  error.   A  discussion  of  the 

statistics  of  prime  interest  is  provided  in  IV. 
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TABLE  1.   Regression  coefficients  for  NESC  dependent  data 

(a)  500  mb ,  three-day  sample 


CONST 

co 

cl 

C2 
C3 
C4 
C5 
C6 
C7 
C8 


C0 
Cl 
C2 
C3 
C4 
C5 
C6 
C7 
C8 


20-40N 

35-55N 

50-70N 

-43.1058 

-52.2372 

-90.9234 

-0.0327 

0.4187 

0.1508 

-0.0932 

-0.4795 

-0.2562 

0.8922 

-0.2670 

0.0 

0.6949 

3.0232 

2.9128 

-1.8871 

-2.6461 

-1.7719 

-0.9188 

-0.0726 

-0.1134 

0.4129 

0.1608 

0.2351 

0.6664 

0.0764 

-0.1635 

(b)  300  mb ,  three-day 

sample 

-70.5523 

-113.6468 

-108.4219 

0.4109 

0.5032 

0.3409 

-0.2953 

-0.3192 

-0.2337 

-1.4766 

-1.4030 

-0.6563 

3.0286 

3.1138 

1.3928 

-0. 7404 

0.0857 

1.0938 

-0.5295 

* 

0.0 

0.1261 

-0.5835 

-1.0623 

-1.1938 

0.6865 

0.4809 

0.3922 

indicates  that  this  variable  (*)  was  not  accepted 
because  it  failed  to  explain  appreciably 
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TABLE  2.   Regression  coefficients  for  NESC  dependent  data 

(a)  500  mb  ,  four-day  sample 


CONST 

20-40N 

35-55N 

50-70N 

co 

-48.6906 

-51.8292 

-92.4865 

cl 

0.0 

0.4842 

0.1143 

C2 

-0.0907 

-0.6454 

-0.1644 

C3 

0.8032 

-0.1472 

-0.0980 

C4 

0.6590 

3.0011 

3.0117 

C5 

-1.5652 

-2.6431 

-1.7844 

C6 

-0.9017 

0.0 

-0.1327 

C7 

0.2990 

0.1605 

0.1894 

cQ 

0.6797 

0.0354 

-0.1311 

(b)  300  mb ,  four-day  sample 
CQ  -101.7840       -108.5859      -108.1667 

ci 

C2 
C3 
C4 
C5 
C6 
C7 
C8 

* 

indicates  that  this  variable  (*)  was  not  accepted 
because  it  failed  to  explain  appreciably 


0.1288 

0.5599 

0.3757 

-0.0794 

-0.5138 

-0.3326 

-0. 2788 

-1.0924 

-0.5185 

0.5219 

2.7352 

1.2617 

1.2700 

0.1432 

1.1646 

-1.0121 

* 
0.0 

0.1026 

-0.5885 

-1.0406 

-1.2316 

1.0295 

0.5160 

0.4386 
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TABLE  3.   Regression  coefficients  for  FNWC  dependent  data 

(a)  500  mb ,  three-day  sample 


CONST 

20-40N 

35-55N 

50-70N 

co 

-15.6121 

-55.3823 

-85.1272 

cl 

0.3605 

0.4119 

9.0638 

C2 

-0.8342 

-0.8714 

-0.7219 

C3 

0.7760 

0.9170 

1.7073 

C4 

1.2602 

1.6287 

0.6259 

S 

-1.8509 

-1.9698 

-0.6344 

C6 

-0.9308 

0.1453 

-0.1520 

C7 

0.2310 

-0.1858 

0.0* 

Co 

0.4440 

0.1834 

-0.0545 

(b)     3.00    mb ,    three-day    sample 


-56.7715 

0.4459 

-0.6807 

-0.3890 

1.9185 

-0.5798 

-0.9211 

-0.2675 

0.5916 


-107.1405 

0.2963 

0.0459 

-1.1144 

2.1166 

0.3357 

-0.1590 

-0.5875 

0.2311 


-99.6836 
0.3651 

-0.6303 
0.3020 
0.5930 
0.9944 
0.0444 

-0.6276 
0.0392 


indicates  that  this  variable  (*)  was  not  accepted 
because  it  failed  to  explain  appreciably 
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TABLE  4.   Regression  coefficients  for  FNWC  dependent  data 

(a)  500  mb  ,  four-day  sample 


CONST 

20-40N 

35-55N 

50-70N 

co 

-15.3487 

-51.3814 

-88.3277 

cl 

0.4477 

0.5272 

0.0260 

C2 

-1.0928 

-1.2317 

-0.5774 

C3 

1.0306 

1.3415 

1.4995 

CA 

1.1163 

1.3383 

0.8176 

C5 

-1.7344 

-1.9340 

-0.6736 

C6 

-0.8981 

0.2096 

-0.1147 

C7 

0.2645 

-0.1328 

-0.0979 

Co 

0.3767 

0.1141 

0.0 

(b)  300  mb ,  four-day  sample 


-53.5409 

0.4850 

-0.  7624 

-0.4110 

2.0702 

-0.6480 

-0.8955 

-0.1607 

0.4292 


-103.0841 

0.3744 

-0.1827 

-0.8859 

2.0795 

0.1929 

-0.1818 

-0.5045 

0.2066 


-99.9952 
0.3821 

-0.6721 
0.3722 
0.5005 
1.0677 
0.0171 

-0.6575 
0.0757 


indicates  that  this  variable  (*)  was  not  accepted 
because  it  failed  to  explain  appreciably 
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IV.   DEPENDENT  DATA  RESULTS 

The  detailed  results  of  the  tests  on  dependent  data  are 
available  in  Tables  5  and  6.   However,  the  essential  por- 
tions needed  for  comparisons  by  the  different  tests  using 
dependent  data  have  been  listed  in  Table  7  where  the  values 
of  fractional  explained  variance  as  originally  listed  in 

Tables  5  and  6  are  summarized.   Also  listed  in  Table  7  are 

2 
the  values  of  R   averaged  over  the  data-base  specification 

2 
period  (or  mode)  for  each  band  as  well  as  the  values  of  R 

averaged  over  the  bands  for  a  given  data-base  mode.   Com- 
parisons between  regressions  were  based  on  the  fractional 

2 
explained  variance  R   rather  than  on  the  values  of  standard 

2 
errors,  as  R   has  been  normalized  with  respect  to  the  sam- 
ple variances,  which  vary  between  bands. 

The  third  value  listed  in  the  band  columns  of  both  parts 
(a)  and  (b)  of  Table  7  is  the  average  over  the  data-base 
modes  in  each  band.   The  values  listed  in  the  mode-average 

row  and  band-average  column  are  the  overall  averages  by  the 

2 
statistical  regression.   Comparisons  of  R   are  made  of  the 

overall  averages  and  of  band  averages,  as  well  as  data-base 

mode  averages  at  both  the  500-and  300-mb  levels.   Differ- 

2 
ences  in  values  of  certain  compared  R   are  termed  "%  shrink- 
age" in  the  explained  variance  arising  from  the  specifica- 
tion tests  under  comparison. 
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A.  COMPARISON  OF  NESC  SPECIFICATION  STATISTICS 
Examination  was  first  made  of  the  left  hand  side  of 

Table  7  for  comparison  of  the  NESC  regression  results  at 

300  mb  relative  to  those  at  500  mb .   The  outstanding 

2 
result  of  this  comparison  reveals  RN(3)  to  be  smaller  than 

2 
R  (5)  by  0.1660,  or  16.6%  when  averaged  over  bands  1,  2, 

and  3.   Upon  comparison  of  individual  band  results  (averaged 

over  the  three  and  four  day  specification  modes)  one  finds 

the  same  order  of  shrinkages  distributed  over  the  bands, 

namely,  an  average  of  16.6%,  but  with  a  larger  than  normal 

2 
shrinkage  in  band  2  where  the  value  of  R„  fell  off  to 

N 

0.6857. 

Upon  use  of  the  four-day  estimation  equation  (10)  in 
the  respective  cases  (left  side  of  Table  7)  there  was  a 
small  shrinkage  in  explained  variance  (relative  to  the 
three-day  results)  at  both  the  300-  and  500-mb  levels. 
This  result  was  not  quite  anticipated  since  Smith   et  aJ 
found  it  advisable  to  use  data-banks  of  two  weeks  for 

specification.   Nevertheless,  the  small  differences  between 

2 
the  R   for  the  two  prediction-base  methods,  whether  at  300 

or  500  mb ,  were  not  large  enough  to  be  considered  signifi- 
cant in  the  tests  conducted  here. 

B.  COMPARISON  OF  FNWC  SPECIFICATION  STATISTICS 

Here  our  comparisons  will  mainly  concern  the  right  hand 
side  of  Table  7.  Again,  predictor-mode  differences  seem  to 
be  marginal  so  that  for  the  present  they  will  be  disregarded 
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2  2 

Comparison  of  Rj.(3)  with  Rp(5)  shows  that  band  average 

2        2 
shrinkage  in  the  amount  R.,(5)  -  R^O)  =  0.1480  occurs.   A 

breakdown  of  this  statistic  by  bands  shows  shrinkage  in  the 
amounts  9.16%,  22.39%,  and  12.87%  for  bands  1,  2,  and  3, 
respectively . 

Comparisons  made  between  NESC  and  FNWC  regression  re- 
sults (left  side  and  right  side  of  Table  7)  show  approxi- 
mately 2%  difference  between  R^CS)  =  0.8737  and  RM(5)  = 

r  N 

0.8957  and  less  than  0.5%  difference  at  the  300-mb  level. 
Thus,  the  FNWC  specifications  at  both  levels  explained  very 
nearly  as  much  variance  as  did  the  NESC  results.   There 
were  some  variations  between  bands,  but  generally,  the 

same  geographic  pattern  persisted  throughout:  namely, 

2 
moderate  to  high  R   in  bands  1  and  3,  with  the  lowest  value 

2 
in  R   occurring  in  band  2.   This  was  true  of  all  FNWC  speci- 
fications (300  and  500  mb )  as  well  as  of  the  NESC  300  mb 

specification.   The  NESC  500  mb  specification  gave  rise  to 

2 
an  increase  in  R   with  an  increase  in  latitude. 

In  summary,  both  NESC  and  FNWC  regression  equations 

tested  well  in  terms  of  explained  variance  at  500  mb .   Both 

showed  diminished  specification  at  300  mb  to  values  of 

2  . 
R  ■  0.73.   Finally,  in  no  case  was  there  a  significant 

change  when  testing  the  four-day  data  base  results  against 

those  for  the  three-day  data  base. 
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TABLE  5.   Dependent-data  temperature  specification 
statistics  resulting  from  pooled  samples  of  NESC  data 


Band 
20-40N 
35-55N 
5 0-7  ON 

20-40N 
35-55N 
50-70N 

20-40N 
35-55N 
50-70N 

20-40N 
35-55N 
50-70N 


(a)  500  mb ,  three-day  sample  base 


Sample 
Size    Mean 


Std. 
Dev. 


171  -15.5029  6.8730 
191  -28.3298  6.9763 
177      -34.1469   7.0397 


R 

0.8533 
0.9113 
0.9321 


(b)  500  mb ,  four-day  sample  base 
206      -15.4369   6.7276    0.8379 
228      -28.1623   6.7376    0.9040 
216      -34.2176   7.0982    0.9358 

(c)  300  mb ,  three-day  sample  base 
171      -42.7895   4.6668    0.7719 
191      -50.4869   4.7427    0.6842 
177      -54.6667   3.9610    0.7833 

(d)  300  mb ,  four-day  sample  base 
206      -42.7233   4.6398    0.6548 
228      -50.3509   4.6875    0.6871 
216      -54.7778   3.9610    0.7970 


Std. 
Error 

2.6970 

2.1218 

1.8724 

2  .2275 
2.1202 
1.8336 

2.2834 
2.7152 
1.8874 

2.7809 
2.6637 
1.8186 


31 


TABLE  6.   Dependent-data  temperature  specification 
statistics  resulting  from  pooled  samples  of  FNWC  data 


Band 
20-40N 
35-55N 
50-70N 

20-40N 
35-55N 
50-70N 

20-40N 
35-55N 
50-70N 

20-40N 
35-55N 
50-70N 


Sample 
Size 


(a)  500  mb ,  three-day  sample  base 

2 


Std. 
Dev. 


R 

0.8760 
0.8704 
0.8816 


Mean 

171  -15.9520   6.6156 

191  -28.2174   7.3512 

177  -34.2333   7.3384 

(b)  500  mb ,  four-day  sample  base 
206  -15.9179   6.4674  0.8651 
228  -28.0130   7.1574  0.8598 
216  -34.2405   7.3470  0.8894 

(c)  300  mb ,  three-day  sample  base 
171  -42.2974   5.2983  0.7833 
191  -50.8692   4.6688  0.6413 
177  -54.9912   4.3548  0.7494 

(d)  300  mb ,  four-day  sample  base 
206  -42.3044   5.1558  0.7747 
228  -50.7638   4.6363  0.6411 
216  -55.1295   4.3142  0.7641 


Std. 
Error 

2.3861 

2.7045 

2.5767 

2.4232 
2  .7283 
2.4842 

2.5267 
2.8571 
2.2313 

2.4962 
2.8276 
2.1354 
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V.   INDEPENDENT  DATA  TEST  RESULTS 

Verification  of  the  dependent  regression  equations  was 
accomplished  using  the  regression  coefficients  (described 
in  III)  and  the  12-and  24-hour  independent  data  samples. 
The  equations  of  form  (10)  developed  from  both  the  three-day 
and  four-day  samples  were  tested.   The  independent  tests 
generate  equations  of  the  general  form 

T(5,3)  =  BQ  +  B1T(5,3)  (13) 

where  Bn  and  B   are  regression  constants  and  T  represents 
the  temperature-estimators  from  the  dependent  regression 
equation  (10)  with  .coefficients  appropriate  to  the  level, 
the  band  and  the  mode-stratifications  (Tables  1,2,3,4). 
The  regression  analysis  for  these  tests  is  similar  to  the 
techniques  applied  in  IV  in  that  T(5,3)  is  compared  by  re- 
gression methods  to  the  observed  temperature  at  each  scan- 
spot.   The  results  of  the  independent  tests  are  given  in 
Tables  8,  9,  10,  and  11.   Table  12  is  constructed  in  a  man- 
ner similar  to  that  used  to  develop  Table  7.   The  major 
difference  between  Tables  12  and  7  is  that  the  12-  and  24- 
hour  independent  samples  are  statistically  pooled  in  appli- 
cation of  the  independent  regression  test  in  order  to  obtain 
a  larger  sample  size  for  each  band,  and  also  for  ease  of 
comparisons.   The  results  of  the  pooled  independent  data 
verifications  are  presented  in  Table  12. 
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Internal  comparisons  of  the  values  of  fractional  ex- 
plained variance  for  the  independent  data  yielded  results 
quite  similar  to  those  obtained  from  the  dependent  data. 
The  comparisons  are  stratified  in  two  manners;  by  bands  and 
by  modes,  as  well  as  by  levels. 

A.   COMPARISONS  OF  NESC  INDEPENDENT  VERIFICATIONS 

From  the  left  side  of  Table  12,  it  is  easily  observed 

2 
that  the  largest  overall  average  R  =  0.8392  occurred  in 

connection  with  the  500-mb  verification  of  T„(5).   The  cor- 

N 

responding  mean  NESC  result  at  300  mb  verified  with  a  per- 

2 
centage  explained  variance  R  (3)  =  0.7533,  i.e.,  with  a 

2 
shrinkage  of  8.59%  relative  to  R  (5). 

As  noted  in  IV(A  and  B),  both  NESC  dependent  specifica- 

2 
tions  (300  and  500  mb)  in  band  3  verified  highest  with  R 

ranging  between  0.82  and  0.91.   The  previously  noted  trend 

2 
in  R   continues  at  the  500-mb  level  of  independent  data  with 

2 
R   increasing  poleward  from  0.78  in  band  1  to  0.91  in  band 

2 
3.   At  300  mb ,  this  band  distribution  of  R  (3)  is  not  so 

apparent  since  values  close  to  0.72  are  associated  with 

2     ~ 
bands  1  and  2  with  a  mean  value  of  R  (3)  ~  0.82  existing  in 

band  3 . 


B.   COMPARISONS  OF  FNWC  INDEPENDENT  VERIFICATIONS 

Comparisons  of  the  FNWC  independent  tests  (Table  12, 
right  side)  bring  out  primarily  the  same  results  which  have 

been  developed  from  dependent  base  samples.   For  example, 

2 
the  FNWC  overall  Rr,(5)  ~  0.78  may  be  compared  with 

r 
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2  2 

R  (3)    0.75.   In  general,  the  independent  R  (3, indep)  is 

2 
close  to  the  values  found  for  R  (3,dep).   In  keeping  with 

i 

earlier  results,  band  2  has  a  relative  minimum  value  of 

2 
R  (indep)  of  magnitude  0.6760.   The  corresponding  results 

2  2 

for  Rp(5, indep)  relative  to  R  (3, indep)  show  that  the  for- 
mer averages,  band-f or-band ,  between  6  and  8%  higher  than 
the  latter.   Comparison  of  FNWC  and  NESC  independent  tests, 
both  at  500  mb ,  shows  higher  mean  specification  of  the 
latter  by  close  to  5%  in  favor  of  NESC  over  FNWC  independent 
verification  in  bands  2  and  3.   The  NESC  verification  sta- 
bility seems  to  be  dominant  in  bands  2  and  3  with  values  of 

2  ~ 
R   ~  0.85.   Here  again,  the  geographic  trend  previously 

noted  is  present. 

2 
As  has  been  described  before,  R  (5, indep)  also  has  a 

r 

2 
band-2  minimum  but  at  a  level  of  10%  higher  than  R  (3, indep) 

r 
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VI.   COMPARISONS  OF  INDEPENDENT  TO  DEPENDENT  RESULTS 

A  summary  of  independent  to  dependent  test  results  is 
given  in  Table  13.   Three  groups  of  tests  are  included. 
The  first  stratification  (part  a,  Table  13)  is  averaged 
over  all  the  bands  and  the  data  base  modes.   The  outstanding 
aspect  of  part  (a),  column  2,  is  that  it  summarizes  the 

comparative  explained  variances  of  T   and  T  .   The  high 

F        N 

explained  variances  for  the  dependent  case  at  500  mb  may  be 

2 
noted,  as  well  as  the  relatively  modest  decreases  in  R 

(between  6  and  8%)  upon  testing  with  the  independent  pooled 

samples . 

Part  (b)  of  Table  13  averages  over  pressure  levels  and 

illustrates  in  summary  form  the  distinction  between  the  NESC 

and  FNWC  regression  samples.   These  samples  show  the  geo- 

2  2 

graphic  patterns  of  R   noted  earlier:  e.g.  R  (dep)  has  a 

2 
tendency  to  increase  poleward  whereas  R  (dep)  has  a  band  2 

r 

minimum.  The  independent  test  column  of  part  (b)  summarizes 

the  pattern  of  T  -verification  by  bands  (but  independent  of 

r 

2 
levels).   Upon  verification,  the  mean-regression  R   shrink- 
age of  T   averages  between  4  and  6%  when  all  bands  are  con- 
F 

s  ider ed . 

Part  (c)  of  Table  13  summarizes  some  of  the  key  ideas 

put  forth  in  earlier  sections.   For  example,  the  variable 

2 
T  (5  dep)  is  characterized  by  a  high  value  of  R  (5).   The 

"latter  statistic  decreases  slightly  (by  about  5%)  when 
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applied  to  independent- tes t  data.   The  same  remarks  may  be 

2 
applied  to  T  (5)  with  slightly  smaller  Rp(5)  ~  0.874  noted 

in  the  dependent  case  coupled  with  a  somewhat  larger  shrink- 
age (about  9%)  after  test  upon  independent  data.   The  con- 
clusions regarding  the  comparative  instability  of  the  depen- 
dent test  for  T  (3)  become  clear  from  Table  13,  part  (c) , 

2 
where  it  is  seen  that  the  values  of  R   increase  approxi- 
mately 3%  from  dependent  to  independent  tests. 

The  results  for  T  (3  dep)  and  T  (3  indep)  are  inconclu- 

F  r 

sive  in  the  form  presented  in  Table  13(c).   The  whole  set 
of  results  of  part  (c)  indicates  that  little  if  any  addi- 
tional specification  accuracy  is  obtained  from  the  four-day 
data  base  in  lieu  of  the  three-day  data  base. 
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TABLE  13.   Summary  of  independent  to  dependent  test  results 


(a) 


(b) 


(c) 


test  being 

dependent  „ 
value  of  R 

compared 

V5) 

TpC5) 

TN(3) 
V3) 

0.8957 
0.8737 

0.7298 

0.7257 

T  ,  band 

Tp,  band 

T*   band 

Tp,  band 

T*   band 

T^,  band 
F 

1 

0.7795 

1 

0.8248 

2 

0.7966 

2 

0.7532 

3 

0.8620 

3 

0.8212 

TN(5) 
TN(5) 
TN(3)" 

Tj(5) 

Tp(5), 

lj(3) 

,  3 

,  4 

3 

day 
day 
day 

0.8989 
0.8928 
0.7465 

4 
3 

day  . 
day 

0. 7129 
0.8760 

4 

day 

0.8714 

,  3 

day 

0.7247 

,  4 

day 

0.7266 

independent 
value  of  R 


0.8392 
0.7798 
0.7533 
0.7201 


0 
0 
0 
0 
0 
0 

0 
0 
0 
0 
0 
0 
0 
0 


7548 
7561 

7701 
7177 
8640 
7761 

8369 
8369 
7841 
7226 

7725 
7871 
7504 
6898 
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VII.   CONCLUSIONS 

The  results  of  Sections  IV,  V,  and  VI  afford  confirming 
evidence.   For  each  case,  high  specification  for  T   occurs 

at  the  500-mb  level.   In  spite  of  considerable  shrinkage, 

2 
there  still  remains  a  significant  specification  value  of  R 

at  300  mb .   This  was  true  also  of  the  FNWC  specifications. 

2 
The  reduction  in  R   associated  with  300  mb  T._  and  T„ 

N       F 

specifications  was  maximal  in  band  2  where  more  inter-di- 
urnal variability  in  the  temperature  and  associated  tropo- 
pause  conditions  occur.   The  NESC  specifications  obtained  by 
Smith  et  al  (1970)  by  the  16-predictor  equation  also  fell 
off  to  a  relative  minimum  at  the  tropopause  level  in  their 
SIRS-A  test. 

Not  only  is  part  of  the  observed  lack  of  specification 
associated  with  the  sampling  of  temperature  extremes  at  or 
near  the  300-mb  level,  but  also  the  inter-diurnal  varia- 
bility in  cloud  cover  would  affect  the  noise  contained  in 
the  radiances  due  to  the  cloud-correction  technique,  which 
is  not  error-free.   The  errors  due  to  cloud-correction  must 
be  maximal  in  channels  2  and  3  for  which  maximum  dl/dln  p 
weights  occur  in  the  upper  troposphere.   Hence,  the  mid- 
latitude  band  (band  2)  which  had  most  inter-diurnal  as 
well  as  geographical  sampling  variability  had  the  poorest 
300-mb  specification.   On  the  other  hand,  the  500-mb  speci- 
fication and  verification  was  adequate  throughout.   The 


45 


fact  that  the  best  300-mb  specification  occurred  in  band  3 
(50-70N)  where  cloud  tops  are  climatologically  lower  than 
in  the  other  two  bands  tends  to  confirm  this  conclusion. 
Similar  tests  at  additional  levels  would  have  provided  a 
better  perspective  in  assessing  the  use  of  the  SIRS  cor- 
rected radiance  data  in  constructing  vertical  profiles.   As 
a  result,  a  more  balanced  view  of  the  vertical  resolution 
afforded  by  the  SIRS  regression  techniques  would  have 
evolved . 

Several  factors  could  work  to  improve  the  operational 
routine  tested  here:  (1)  to  make  use  of  a  multi-level 
cloud-correction  model  for  correction  of  SIRS-radiances 
(similar  to  that  described  by  Smith  et  al    1970);  (2)  to 
employ  a  16-predictor  equation  involving  squares  of  radi- 
ances as  well  as  linear  terms;  (3)  to  use  a  larger  sample 
base  confined  to  a  smaller  band  zone,  for  example,  15 
latitude  overlapping  bands  rather  than  the  20   bands  used 
here  . 

It  is  felt  that  with  these  modifications,  a  SIRS 
operational  unit  attached  to  FNWC  could  profitably  use  the 
SIRS  data  in  improving  their  update  analysis  over  sparse 
data  regions,  and  in  determining  the  onset  of  significant 
circulation  changes. 
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